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1 Abstract 

Why do some clusters have cool cores while others do not? In this paper, 



> 

o 

\ cosmological simulations, including radiative cooling and heating, are used 

■ to examine the formation and evolution of cool core (CC) and non-cool core 
(NCC) clusters. Numerical CC clusters at z — accreted mass more slowly 
over time and grew enhanced cool cores via hierarchical mergers; when late 
major mergers occurred, the CCs survived the collisions. By contrast, NCC 

■ clusters of similar mass experienced major mergers early in their evolution 
Onj that destroyed embryonic cool cores and produced conditions that prevent 
q ■ CC re-formation. We discuss observational consequences. 

2 Introduction 

k> , Recently, 49% of clusters were identified as having cool cores in a flux-limited 

sample, HIFLUGCS, based upon ROSAT and ASCA observations Q. The 
earliest and simplest model assumed these clusters to be spherical, isolated 
systems where "cooling flows" formed; as radiating gas loses pressure support, 
cooling gas flows inwards to higher density values which further accelerates 
the cooling rate (e.g., |2|). However, the predicted end-products of this mass 
infall (e.g., star formation, HI, CO) have not been observed and the central 
temperatures indicate that the gas at the cores has cooled only moderately 
(generally only 30-40% of the virial temperature |3j). Furthermore, this model 
failed to consider the important effects of mergers and on-going mass accretion 
from the surrounding supercluster environment. 
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We have performed numerical simulations of the formation and evolution 
of clusters in a cosmological context using the adaptive mesh refinement N- 
body/hydro code Enzo, aimed at further understanding cool cores Q- We 
find that most mergers are oblique with halos spiraling into the cluster cen- 
ter and gently bequeathing cool gas; this allows the cool core to strengthen 
over time. Thus, in this model, cool cores themselves grow hierarchically via 
the merger/accretion process. Also, a realistic model of cool cores must in- 
volve heating, as well as cooling, that "softens" the cores (i.e., reduces density 
and increases temperature) potentially by star formation or by AGNs. This 
softening would increase the susceptibility of disruption of the cool core by 
subsequent mergers. 

In this paper, we present new enhanced-physics Enzo cosmology simula- 
tions that include radiative cooling, star formation (i.e., a mass sink for cold 
gas), and heating (modeled via kinetic energy input from Type II supernovae). 
In addition to the more realistic baryonic physics, these simulations are supe- 
rior to previous numerical models (e.g., |i|(|) due to their bigger volumes and 
larger samples of clusters. With hundreds of clusters of mass > 10 14 M©, we 
have the capability for the first time to examine evolutionary effects in these 
numerical clusters with adequate statistics. 



3 The Formation of CC and NCC Clusters 

Enzo is an adaptive mesh refinement cosmological simulation code that cou- 
ples an N-body particle-mesh (PM) solver with an Eulerian block-structured 
AMR method for ideal gas dynamics. Enzo has physics modules which model 
radiative cooling of both primordial and metal-enriched gases, and prescrip- 
tions for star formation and feedback [tJ. 

Our AMR simulations have produced a sample of 1494 numerical clusters 
with < z < 1 and M V i r i a i > 10 14 Af Q . Only «10% of this sample of numerical 
clusters have cool cores. This fraction is smaller than observed. We believe 
that our specific recipe for star formation and heating has over-softened the 
cores making them too vulnerable to disruption. This will be addressed in 
subsequent simulations. 

In Figure 1, we show examples of the evolution of a CC and an NCC 
cluster which represent well the general scenarios for how such clusters form. 
This figure illustrates the different histories of two clusters which have the 
same final mass at z = 0, M V i r i a i — 5 x 10 14 M Q . 

Many clusters start off with cool cores early in their history when ini- 
tial conditions produce central densities and temperatures that allow the gas 
to radiatively cool. Depending upon the sequence and magnitude of cluster 
collisions, cool cores will be destroyed in some clusters and preserved in oth- 
ers. NCC clusters often undergo major mergers early in their history. From 
z=0.65 to z=0.5, the mass of the NCC cluster in Fig. 1 increased by 110%. 
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Fig. 1. The time history of the evolution of a CC (left) and an NCC (right) numerical 
cluster. Dashed lines indicate the total mass in units of the virial mass normalized to 
1 at z—0. Solid line is the central temperature normalized by the virial temperature. 
Note that the dip in mass at z ~ 0.15 in the NCC cluster is an artifact of how 
we measure the cluster mass within the virial radius; it is caused by the passage 
of an infalling subcluster through the core and back out beyond r V i T i a i with a final 
acquisition into the cluster by z—Q. 



Such mergers destroy the cool cores, leaving behind hotter, thermalized, mod- 
erately dense cores. Subsequently, cool halos infalling into these NCC clusters 
do not survive passage through the central parts of the clusters and the cen- 
tral conditions do not allow cool cores to re-establish. NCC clusters continue 
to experience minor mergers as they continue to slowly evolve (typically mass 
increases 10-20% over Gyr time frames). We suggest that such an early major 
merger produced the complex characteristics observed today for the Coma 
cluster. 

On the other hand, Figure 1 suggests that CC clusters evolve differently. 
Many CC clusters do not experience major mergers early in their history but 
rather grow slowly such that the cool cores increase in mass and stability. 
The CC cluster at z=0.65 began with a higher mass than the NCC at the 
same redshift. From z=0.65 to z=0.1, mergers increased the mass of the CC 
cluster by only «20%. At z ss 0.1, substantial multiple mergers occurred 
increasing the mass by ~50% by z=0. But by this epoch, the cool core is so 
well established that it survives ever larger impacts. In fact, the cool gas from 
the infalling halos near the present epoch merges with the cool core such that 
it continues to grow. 



4 Consequences of Evolutionary Differences in CC and 
NCC Clusters 

The large-scale X-ray surface brightness profiles (Sx) of our simulated NCC 
clusters have Sx distributions which fit very well to /3-models (Sx oc [1 + 
(r/r c ) 2 ] 1 / 2 ~ 3 P) out to the cluster virial radius as we recently described 0. 
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This is not true for CC clusters in our samples. As detailed in |8(, the (3- 
model does a poor job of extrapolating the surface brightness from r 50 o to 
^200 (subscripts refer to overdensities relative to critical density; r 2 oo ~ r viria i) 
for CC clusters. The CC cluster profile has a typically «35% steeper slope 
than the /3-model in the outer part of the cluster. We view this as a further 
signature of recent merger activity. The result is a bias or overestimate of the 
density profile, leading to a significant overestimate in mass. 
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Fig. 2. Histograms of the quality of fit, as measured by the reduced \ 2 , f° r /3-model 
fits to the X-ray surface brightness profiles for ~800 of our numerical clusters with 
M > 5 x 10 14 M Q . Dotted line is NCC clusters and solid line is CC clusters. 



These results are typical for CC clusters in our sample. We have fit (3- 
models to the synthetic X-ray surface brightness profiles for all the numerical 
clusters with < z < 0.5 and M > 5 x 10 14 M Q («800 clusters). The fits 
excluded the core regions with r < 100 hr 1 kpc and extended out to rsoo to 
mimic what is typically available for X-ray observations. We then extrapolated 
the fit out to r2oo and calculated the reduced \ 2 between the /3-model and 
the actual profile for the numerical clusters. A histogram of the cumlative 
fraction of the \ 2 distribution for NCC and CC clusters is presented in Figure 
2. There is a significant difference between the quality of the /3-model fits 
between these two types of clusters. Nearly 90% of the NCC clusters have 
reduced \ 2 < 1 whereas only 60% of CC clusters have values < 1. This results 
in a mass overestimate or bias of typically a factor of 3-5 for CC clusters. Thus, 
there is a breakdown in the assumptions of a single /3-model (i.e., dynamical 
equilibrium as assumed in simple cooling flow model) even outside the core 
for CC clusters. Similar results for a sample of cool core clusters observed by 
Chandra have been found || . 

Our numerical CC clusters also generally show a higher incidence of cold 
substructure outside the core than do their NCC counterparts when measured 
relative to the virial temperatures. The lower temperature gas is composed 
of both compact infalling halos as well as diffuse cool gas. We predict that 
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this signature will be apparent in hardness ratio maps. As shown in Fig. 3 
for two typical cases drawn from our simulations, the hard-to-soft band ratios 
(2-8 keV/0.7-2 keV) illustrate the abundance of cooler gas beyond the cores in 
CC clusters. We find that this prediction agrees well with the hardness ratio 
distribution of several Abell clusters observed by Chandra. 




) 2-8kev/So.5-2keV 



Fig. 3. Maps (3 h' 1 Mpc on a side) of the X-ray hardness ratio (2-8 keV)/(0.7-2 
keV) for numerical clusters with a cool core (left) and without a cool core (right). 
Dark are smaller (0.2) HR values (corresponding to cool temperatures) and white 
are large (> 0.8) HR values (hot temperatures). The clusters have similar masses of 
?s 5 x 10 14 Mq . At the bottom, the histogram is the distribution of hardness ratios for 
the two clusters above; dotted is for the NCC cluster and solid is for the CC cluster. 
The central 100 kpc of each cluster (cool core region) was excluded in constructing 
this histogram. 



5 Conclusions 

Our cosmological numerical simulations suggest that 
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• Non-cool core (NCC) clusters suffer early major mergers when embryonic 
cool cores are destroyed. Cool core (CC) clusters, on the other hand, grow 
more slowly without early major mergers. 

• The X-ray surface brightness profiles for NCC clusters are well fit by sin- 
gle /3-models whereas the outer emission for CC clusters is biased low 
compared to /3-models. The resulting densities and masses of CC clusters 
estimated from single /3-model extrapolations are biased high by factors 
of 3-5. 

• CC clusters have more cool gas beyond the cores than do NCC clusters, 
reflected in hardness ratio maps. 

This work was supported in part by grants from the U.S. National Science 
Foundation (AST-0407368) and NASA (TM3-4008A). 
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